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Preface 


This text contains a comprehensive discussion on continuous and discrete time signals and 
systems with many MATLAB® and several Simulink® examples. It is written for junior and 
senior electrical and computer engineering students, and for self-study by working professionals. 
The prerequisites are a basic course in differential and integral calculus, and basic electric circuit 
theory. 

This book can be used in a two-quarter, or one semester course. This author has taught the 
subject material for many years and was able to cover all material in 16 weeks, with 2 Vi lecture 
hours per week. 

To get the most out of this text, it is highly recommended that Appendix A is thoroughly 
reviewed. This appendix serves as an introduction to MATLAB, and is intended for those who 
are not familiar with it. The Student Edition of MATLAB is an inexpensive, and yet a very 
powerful software package; it can be found in many college bookstores, or can be obtained directly 
from 

The MathWorks™ Inc., 3 Apple Hill Drive, Natick, MA 01760-2098 
Phone: 508 647-7000, Fax: 508 647-7001 
http://www .mathworks. com 
e-mail: info@mathworks.com 

The elementary signals are reviewed in Chapter 1, and several examples are given. The purpose of 
this chapter is to enable the reader to express any waveform in terms of the unit step function, and 
subsequently the derivation of the Laplace transform of it. Chapters 2 through 4 are devoted to 
Laplace transformation and circuit analysis using this transform. Chapter 5 is an introduction to 
state-space and contains many illustrative examples. Chapter 6 discusses the impulse response. 
Chapters 7 and 8 are devoted to Lourier series and transform respectively. Chapter 9 introduces 
discrete-time signals and the 2 transform. Considerable time was spent on Chapter 10 to present 
the Discrete Lourier transform and LET with the simplest possible explanations. Chapter 11 
contains a thorough discussion to analog and digital filters analysis and design procedures. As 
mentioned above, Appendix A is an introduction to MATLAB. Appendix B is an introduction to 
Simulink, Appendix C contains a review of complex numbers, and Appendix D is an introduction 
to matrix theory. 

New to the Second Edition 

This is an extensive revision of the first edition. The most notable change is the inclusion of the 
solutions to all exercises at the end of each chapter. It is in response to many readers who 
expressed a desire to obtain the solutions in order to check their solutions to those of the author 
and thereby enhancing their knowledge. Another reason is that this text is written also for self- 




study by practicing engineers who need a review before taking more advanced courses such as 
digital image processing. 

Another major change is the addition of a rather comprehensive summary at the end of each 
chapter. Hopefully, this will be a valuable aid to instructors for preparation of view foils for 
presenting the material to their class. 

New to the Third Edition 

The most notable change is the inclusion of Simulink modeling examples. The pages where they 
appear can be found in the Table of Contents section of this text. Another change is the 
improvement of the plots generated by the latest revisions of the MATLAB® Student Version, 
Release 14. 

The author wishes to express his gratitude to the staff of The MathWorks™, the developers of 
MATLAB® and Simulink®, especially to Ms. Courtney Esposito, for the encouragement and 
unlimited support they have provided me with during the production of this text. 

Our heartfelt thanks also to Ms. Sally Wright, RE., of Renewable Energy Research Laboratory 
University of Massachusetts, Amherst, for bringing some errors on the previous editions to our 
attention. 

New to the Fourth Edition 

The most notable change is the inclusion of Appendix E on window functions. The plots were 
generated generated with the latest revisions of the MATLAB® R2008a edition. Also, two end- 
of- chapter exercises were added in Chapter 10 to illustrate the use of the fft and ifft MATLAB 
functions 

The author wishes to express his gratitude to the staff of The MathWorks™, the developers of 
MATLAB® and Simulink®, especially to The MathWorks™ Book Program Team, for the 
encouragement and unlimited support they have provided me with during the production of this 
and all other texts by this publisher. 


Orchard Publications 

www.orchardpublications.com 

info@orchardpublications.com 
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Chapter 1 Elementary Signals 


We can express (1.3) by the waveform shown in Figure 1.2. 


v s 

0 


v 


out 


> t 


Figure 1.2. Waveform for v out as defined in relation (1.3) 


The waveform of Figure 1.2 is an example of a discontinuous function. A function is said to be dis¬ 
continuous if it exhibits points of discontinuity, that is, the function jumps from one value to 
another without taking on any intermediate values. 


1.2 The Unit Step Function u 0 (t) 

A well known discontinuous function is the unit step function u 0 (t) which is defined as 


[0 

t < 0 

u o(t) = j ! 

t > 0 


It is also represented by the waveform of Figure 1.3. 


| Up(t) 


Figure 1.3. Waveform for u 0 (t) 


(1.4) 


In the waveform of Figure 1.3, the unit step function u 0 (t) changes abruptly from 0 to 1 at 
t = 0. But if it changes at t = t 0 instead, it is denoted as u 0 (t -1 0 ). In this case, its waveform and 
definition are as shown in Figure 1.4 and relation (1.5) respectively. 


i — 


u 0 (t-t 0 ) 


0 


Figure 1.4. Waveform for u 0 (t -1 0 ) 


* 


In some books, the unit step function is denoted as u(t) , that is, without the subscript 0. In this text, however, we 
will reserve the u(t) designation for any input when we will discuss state variables in Chapter 5. 
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The Unit Step Function 


Thus, the pulse of Figure 1.9(a) is the sum of the unit step functions of Figures 1.9(b) and 1.9(c) 
and it is represented as u 0 (t) - u 0 (t - 1). 

The unit step function offers a convenient method of describing the sudden application of a volt¬ 
age or current source. For example, a constant voltage source of 24 V applied at t = 0, can be 
denoted as 24u 0 (t) V. Likewise, a sinusoidal voltage source v(t) = V m coscot V that is applied to 
a circuit at t = t 0 , can be described as v(t) = (V m coscot)u 0 (t-t 0 ) V. Also, if the excitation in a 

circuit is a rectangular, or triangular, or sawtooth, or any other recurring pulse, it can be repre¬ 
sented as a sum (difference) of unit step functions. 


Example 1.2 

Express the square waveform of Figure 1.10 as a sum of unit step functions. The vertical dotted 
lines indicate the discontinuities at T, 2T, 3T, and so on. 


A 

v(t) 


® I (3) 

|T \ 2T 

3T 

-A 

0 ! ! 

1 ® 1 

© 




Figure 1.10. Square waveform for Example 1.2 


Solution: 

Line segment ® has height A, starts at t = 0, and terminates at t = T. Then, as in Example 1.1, this 
segment is expressed as 

Vl (t) = A[u 0 (t)-u 0 (t-T)] (1.8) 

Line segment ® has height -A, starts at t = T and terminates at t = 2T. This segment is 
expressed as 

v 2 (t) = -A[u 0 (t-T)-u 0 (t-2T)] (1.9) 

Line segment ® has height A, starts at t = 2T and terminates at t = 3T. This segment is expressed 
as 

v 3 (t) - A[u 0 (t-2T)-u 0 (t-3T)] (1.10) 

Line segment ® has height -A, starts at t = 3T, and terminates at t = 4T. It is expressed as 

v 4 (t) = -A[u 0 (t-3T)-u 0 (t-4T)] (1.11) 
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The Delta Function 


J f(x)8(x)dx - J f(0)5(x)dx + J [f(x) - f(0)]8(x)dx 


(1.39) 


The first integral on the right side of (1.39) contains the constant term f(0); this can be written 
outside the integral, that is, 

j' f(0)8(x)dx = f(0)J l 5(x)dx (1.40) 

—OO —oo 

The second integral of the right side of (1.39) is always zero because 

8(t) = 0 for t<0 and t>0 


and 


Therefore, (1.39) reduces to 


[f(x)-f(0)]| T = 0 = f(0) - f(0) = 0 


J f(x)8(x)dx = f(0) J 8(x)dx 

—oo —oo 

Differentiating both sides of (1.41), and replacing x with t, we obtain 


f(t)5(t) = f(0)S(t) 
Sampling Property of 8(t) 


(1.41) 


(1.42) 


1.4.2 The Sifting Property of the Delta Function 8(t) 

The sifting property of the delta function states that 


J f(t)8(t — ot)dt = f(a) 


(1.43) 


that is, if we multiply any function f(t) by 8(t - a), and integrate from to +°° , we will obtain 
the value of f(t) evaluated at t = a. 

Proof: 

Let us consider the integral 

r b 

I f(t)8(t-a)dt where a<a<b (1-44) 


We will use integration by parts to evaluate this integral. We recall from the derivative of prod¬ 
ucts that 

d(xy) = xdy + ydx or xdy = d(xy)-ydx (1-45) 

and integrating both sides we obtain 
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Chapter 1 Elementary Signals 


6. To drag a line segment to a new position, we place the mouse cursor over that line segment and 
the cursor shape shows the position in which we can drag the segment. 

7. To drag a point along the y-axis, we move the mouse cursor over that point, and the cursor 
changes to a circle indicating that we can drag that point. Then, we can move that point in a 
direction parallel to the x-axis. 

8. To drag a point along the x-axis, we select that point, and we hold down the Shift key while 
dragging that point. 

9. When we select a line segment on the time axis (x-axis) we observe that at the lower end of 
the waveform display window the Left Point and Right Point fields become visible. We can 
then reshape the given waveform by specifying the Time (T) and Amplitude (Y) points. 



Figure 1.26. Waveforms for the Simulink model of Figure 1.24 

The two positive spikes that occur at t = 2, and t = 7, are clearly shown in Figure 1.26. 
MATLAB* has built -in functions for the unit step, and the delta functions. These are denoted by 
the names of the mathematicians who used them in their work. The unit step function u 0 (t) is 

referred to as Heaviside(t), and the delta function 8(t) is referred to as Dirac(t). Their use is illus¬ 
trated with the examples below. 

syms k a t; % Define symbolic variables 

u=k*sym('Fleaviside(t-a) 1 ) % Create unit step function at t = a 

u = 

k*Heaviside(t-a) 

d=diff(u) % Compute the derivative of the unit step function 

d = 

k*Dirac(t-a) 


* An introduction to MATLAB* is given in Appendix A. 
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Chapter 2 


The Laplace Transformation 


T his chapter begins with an introduction to the Laplace transformation, definitions, and 
properties of the Laplace transformation. The initial value and final value theorems are also 
discussed and proved. It continues with the derivation of the Laplace transform of common 
functions of time, and concludes with the derivation of the Laplace transforms of common wave¬ 
forms. 


2.1 Definition of the Laplace Transformation 

The two-sided or bilateral Laplace Transform pair is defined as 


oo 

Sg{f(t)}= F(s) = J f(t)e" st dt 

—oo 


^-‘{F(s)>= f(.) = M 


a + jgd 


a-jo 


F(s)e st ds 


( 2 . 1 ) 


( 2 . 2 ) 


where Sg{f(t)} denotes the Laplace transform of the time function f(t), Sg *{F(s)} denotes the 
Inverse Laplace transform, and s is a complex variable whose real part is a, and imaginary part 
oo, that is, s = a + joo. 

In most problems, we are concerned with values of time t greater than some reference time, say 
t = t 0 = 0, and since the initial conditions are generally known, the two-sided Laplace trans¬ 
form pair of (2.1) and (2.2) simplifies to the unilateral or one-sided Laplace transform defined as 


Sg{f(t)}= F(s) = [ f(t)e st dt - [ f(t)e s 'dt 

\ J o 


(2.3) 


(2.4) 


The Laplace Transform of (2.3) has meaning only if the integral converges (reaches a limit), that 
is, if 
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Chapter 2 The Laplace Transformation 


2.2.3 Frequency Shifting Property 

The frequency shifting property states that if we multiply a time domain function f(t) by an expo¬ 
nential function e at where a is an arbitrary positive constant, this multiplication will produce a 
shift of the s variable in the complex frequency domain by a units. Thus, 


e at f(t) <=> F(s + a) 

Proof: 


(2.14) 


oo oo 

Sg (e at f(t)} = f e at f(t)e s *dt = f f(t)e (s + a)t dt = F(s + a) 

J o J o 


Note 2: 

A change of scale is represented by multiplication of the time variable t by a positive scaling fac¬ 
tor a. Thus, the function f(t) after scaling the time axis, becomes f(at). 


2.2.4 Scaling Property 

Let a be an arbitrary positive constant; then, the scaling property states that 


Proof: 


and letting t = x/a, we obtain 




Sg{f(at)} = [ f(at)e st dt 

J o 


(2.15) 


Sg{f(at)} = f f(x)e~ s(T/a) df- ) = -f f(T)e‘ (s/a)T d(x) = -F ( - ) 

Jq v a ^ a J q a v a / 

Note 3: 

Generally, the initial value of f(t) is taken at t = 0 to include any discontinuity that may be 
present at t = 0. If it is known that no such discontinuity exists at t = 0~, we simply interpret 
f(0") as f(0). 


2.2.5 Differentiation in Time Domain Property 

The differentiation in time domain property states that differentiation in the time domain corre¬ 
sponds to multiplication by s in the complex frequency domain, minus the initial value of f(t) at 

t = 0". Thus, 
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Properties and Theorems of the Laplace Transform 


Proof: 

The Laplace transform 

S£{f(t)} - 


J T f(t)e“ st dt 

f ( t + nT )^ ° -sT 

1 - e 


of a periodic function can be expressed as 

T .2T _3T 

V 


00 1 2 1 

[ f(t)e st dt = [ f(t)e st dt+ [ f(t)e st dt+ [ f(t)e st dt + 
J Q J q J t J 2t 


(2.28) 


In the first integral of the right side, we let t = x, in the second t = x + T, in the third 
t = x + 2T, and so on. The areas under each period of f(t) are equal, and thus the upper and 
lower limits of integration are the same for each integral. Then, 

T T T 

Sg{f(t)} = f f(x)e -ST dx + f f(x + T)e~ s(T + T) dx + f f(x + 2T)e" s(x + 2T) dx + ... (2.29) 

J o J o J o 


Since the function is periodic, i.e., f(x) = f(x + T) = f(x + 2T) = ... = f(x + nT), we can write 
(2.29) as 

T 

S2 {f(x)} - (1 +e~ sT + e~ 2sT + ...)[ f(x)e~ ST dx (2.30) 

J o 

By application of the binomial theorem, that is, 

1 + a + a" + a 3 + ... = —i— (2.31) 

1 - a 

we find that expression (2.30) reduces to 

| f(x)e ST dx 

^ {f(T) } = 

1 - e 


2.2.10 Initial Value Theorem 


The initial value theorem states that the initial value f(0 ) of the time function f(t) can be found 
from its Laplace transform multiplied by s and letting s —> °° .That is, 


limf(t) = lim sF(s) = f(0 ) 

t->o s ^°° 


(2.32) 
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Chapter 2 The Laplace Transformation 


8. Derive the Laplace transform for the sawtooth waveform fsT (t) below. 



9. Derive the Laplace transform for the full-rectified waveform f FR (t) below. 



Write a simple MATLAB script that will produce the waveform above. 
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Partial Fraction Expansion 


F!(s) = 


3s + 2 
s 2 + 3s + 2 


Solution: 

Using (3.6), we obtain 


F^s) = 


3s + 2 


l i 


+ ■ 


The residues are 


and 


3s + 2 = 

s 2 + 3 s + 2 (s + 1)(s + 2) (s + 1) (s + 2) 
3s + 2 


r, = lim (s+ l)F(s) = 

s — > -l (s + 2) 


r 2 = lim (s + 2)F(s) = 3s + 2 
s —> -2 (s+1) 


= -1 


s = -1 


= 4 


s = -2 


Therefore, we express (3.9) as 


F^s) = 


3s + 2 


- 1 + 4 


s+3s + 2 ( s+1 ) (s + 2) 

and from Table 2.2, Chapter 2, Page 2-22, we find that 

1 


-at ., \ 

e u 0 (t) <=> 


s + a 


Therefore, 


F,(s) ---- 1 - —-— <=> (- e ' + 4e 2t )u 0 (t) - f,(t) 

1 (s+ 1) (s + 2) v ’ ow 1W 


(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


The residues and poles of a rational function of polynomials such as (3.8), can be found easily 
using the MATLAB residue(a,b) function. For this example, we use the script 

Ns = [3, 2]; Ds = [1,3, 2]; [r, p, k] = residue(Ns, Ds) 

and MATLAB returns the values 
r = 

4 

-1 

P = 

-2 

-1 

k = 

[] 
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Chapter 3 The Inverse Laplace Transformation 


F 7 (s) - 

Solution: 

Using Steps 1 through 3 above, we obtain 

-2s+ 4 


■ 2s + 4 


F 7 (s) = 


(S 2 +1)(S-1) 2 


rjS + A r 


+ 


21 


+ • 


l 22 


(S 2 +I)(s-1) 2 (S 2 +1) (s-1) 2 


With Step 4, 


■2s + 4 = (r t s + A)(s - l) 2 + r 21 (s“ + 1) + r 22 (s - l)(s“ + 1) 


(3.56) 


(3.57) 

(3.58) 


and with Step 5, 


3 2 

- 2s + 4 = (rj + r 22 )s~ + (- 2r l + A - r 22 + r 21 )s 
+ (rj - 2A + r 22 ) s + (A - r 22 + r 21 ) 


(3.59) 


Relation (3.59) will be an identity is s if each power of s is the same on both sides of this relation. 
Therefore, we equate like powers of s and we obtain 


0 = lq + r 2 , 

0 = -2rj + A-r 22 + r 21 
-2 = r ! - 2A + r 22 
4 = A - r 22 + r 21 

Subtracting the second equation of (3.60) from the fourth, we obtain 

4 = 2r 1 
or 

r i = 2 

By substitution of (3.61) into the first equation of (3.60), we obtain 

0 — 2 + ^22 
or 

r 22 = ~2 

Next, substitution of (3.61) and (3.62) into the third equation of (3.60) yields 


(3.60) 


(3.61) 


(3.62) 


-2 = 2-2A-2 
or 

A = 1 (3.63) 

Finally by substitution of (3.61), (3.62), and (3.63) into the fourth equation of (3.60), we obtain 
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Chapter 3 The Inverse Laplace Transformation 


4/3 , 2/3 2 -4t -t/4, 

-+-<=> r(2e + e ) 

s + 4 s+1/4 3 V 

2 . , 2 


s + 8s + 24s+ 32 (s + 4)(s~ + 4s + 8) (s +4s + 8) ,, , ,• • • 

-= ^-—-- = i-- and by long division 


d. 


s + 6s + ' 


(s + 2)(s + 4) 


s +4s + 8 = s + 2 + — 8'(t) + 28(t) + 4e 

s + 2 s + 2 


(s + 2) 

■2t 


e. 


-2s 3 


F(s)= 


(2s+ 3) 
3 


e 2s F(s)^f(t-2)u 0 (t-2) 

,3 


3/2 


3/8 


3/8 


(2s+ 3) 3 (2 s + 3) 3 /2 3 [(2s + 3)/2] 3 (s + 3/2) 3 


3/12 -(3/2)t\ 3 ,2 -(3/2)t 

«i2! ,e J = l6* e 


e~ 2s F(s)= e- 2 ‘—i-,«^(t-2) 2 e- (3/2Mt - 2, u 0 (t-2) 


(2s + 3) 2 16 


4. The initial value theorem states that lim f(t) = lim sF(s). Then, 

t —» 0 s —> 


f(0) - lim s—— + ^ - - lim + 3s 


s_>00 s + 4.25s + 1 
— lim 


2s"/s" + 3s/s 


s -*“ S + 4.25s + 1 
2 


= lim 


2 + 3/s 


s 2 /s 2 + 4.25s/s 2 + 1/s 2 1 + 4.25/s + 1/s 2 


= 2 


The value f(0) = 2 is the same as in the time domain expression that we found in Exercise 
3(c). 

5. We are given that F(s) = ~ ^ and lim f(t) = lim sF(s) = 10. Then, 

s(s + 1) t—>«> s —> o 

lim s^r^—rr = A lim ——— - -A = 10 


Therefore, 


s— >0 s(s + 1 ) s —> 0 (s + 1) 

F(s) = ~ 10(s ~ 1) = !l + JL = 10_JO^«(io_ 20e" l )u 0 (t) 

s(s + 1) s s + 1 s s + 1 


that is, 

and we observe that 


f(t) = (10 - 20e~ t )u 0 (t) 


lim f(t) - 10 

t —> oo 
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Chapter 4 


Circuit Analysis with Laplace Transforms 


T his chapter presents applications of the Laplace transform. Several examples are presented 
to illustrate how the Laplace transformation is applied to circuit analysis. Complex imped¬ 
ance, complex admittance, and transfer functions are also defined. 


4.1 Circuit Transformation from Time to Complex Frequency 

In this section we will show the voltage-current relationships for the three elementary circuit 
networks, i.e., resistive, inductive, and capacitive in the time and complex frequency domains. 
They are shown in Subsections 4.1.1 through 4.1.3 below. 


4.1.1 Resistive Network Transformation 

The time and complex frequency domains for purely resistive networks are shown in Figure 4.1. 


Time Domain 


+ 


V R (t) = Ri R(t) 

v R (t) R|i R (t) . v R (t) 
■r (,) = — 


+ 


Complex Frequency Domain 

V R (s) - RI r (s) 


v 


V R (s) fC> !r( s ) Vr(s) 

r(S) “ R 


Figure 4.1. Resistive network in time domain and complex frequency domain 


4.1.2 Inductive Network Transformation 

The time and complex frequency domains for purely inductive networks are shown in Figure 4.2. 


Time Domain Complex Frequency Domain 

+ 

v L (t) 


Figure 4.2. Inductive network in time domain and complex frequency domain 

4.1.3 Capacitive Network Transformation 

The time and complex frequency domains for purely capacitive networks are shown in Figure 4.3. 


1 


i L (t) v L (t) = L“ 

iao = l f 


+ 


V L (s) 


4 


sL a It (s) 


V L (s) = sLI l (s) - Li L (0 ) 


v L dt 


Q Li L (0 ) 


I L (s) = 


V L ( S ) , i L (0") 

Ls s 
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Using the Simulink Transfer Fen Block 



Figure 4.23. Plot of v out (t) for Example 4.8. 



Creator 

Figure 4.24. Simulink model for Example 4.8 



Figure 4.25. Waveform for the Simulink model of Figure 4.24 
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Solutions to End-of-Chapter Exercises 


Magnitude Vout/Vin vs. Radian Frequency 



The plot above indicates that this circuit is a second-order low-pass filter. 
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Chapter 5 


State Variables and State Equations 


T his chapter is an introduction to state variables and state equations as they apply in circuit 
analysis. The state transition matrix is defined, and the state-space to transfer function 
equivalence is presented. Several examples are presented to illustrate their application. 

5.1 Expressing Differential Equations in State Equation Form 

As we know, when we apply Kirchoff s Current Law (KCL) or Kirchoff s Voltage Law (KVL) in 
networks that contain energy-storing devices, we obtain integro-differential equations. Also, 
when a network contains just one such device (capacitor or inductor), it is said to be a first-order 
circuit. If it contains two such devices, it is said to be second-order circuit, and so on. Thus, a first 
order linear, time-invariant circuit can be described by a differential equation of the form 

a 'S + a ° y ( t ) = X ( t ) C 5 - 1 ) 

A second order circuit can be described by a second-order differential equation of the same form 
as (5.1) where the highest order is a second derivative. 

An nth-order differential equation can be resolved to n first-order simultaneous differential 
equations with a set of auxiliary variables called state variables. The resulting first-order differen¬ 
tial equations are called state-space equations, or simply state equations. These equations can be 
obtained either from the nth-order differential equation, or directly from the network, provided 
that the state variables are chosen appropriately. The state variable method offers the advantage 
that it can also be used with non-linear and time-varying devices. However, our discussion will 
be limited to linear, time-invariant circuits. 

State equations can also be solved with numerical methods such as Taylor series and Runge- 

Kutta methods, but these will not be discussed in this text . The state variable method is best 
illustrated with several examples presented in this chapter. 


Example 5.1 


A series RLC circuit with excitation 


v s (t) = e 


jcot 


(5.2) 


* These are discussed in “Numerical Analysis Using MATLAB and Excel", Third Edition, ISBN 978-1-934404-03-4- 
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Chapter 5 State Variables and State Equations 



Solution: 


Figure 5.4. Circuit for Example 5.4 


This circuit contains only one energy-storing device, the inductor; therefore, we need only one 
state variable. We choose the state variable to denote the current through the inductor as shown 
in Figure 5.5. 


R 

-WV- 


©T 



(t) = x Ja l 


v s u 0 (0 


By KVL, 


or 


or 


Figure 5.5. Circuit for Example 5.4 with assigned state variable x 
v R + v l = v s u 0 (t) 


Ri + L ft = V s u o« 


Therefore, the state equations are 


Rx + Lx = v s u 0 (t) 


R 1 

X = --x + -v s u 0 (t) 

y = x 


(5.21) 


5.2 Solution of Single State Equations 

If a circuit contains only one energy-storing device, the state equations are written as 

x = ax + (3u 
y = kjX + k 2 u 


(5.22) 


where a, (3, k x , and k 2 are scalar constants, and the initial condition, if non-zero, is denoted as 


x 0 = x(t 0 ) 


(5.23) 
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Chapter 5 State Variables and State Equations 


J 2 2 + l 2 + 3 2 = JlA 
a/i 2 + r + 2 2 = 76 

V(- 1) 2 + 1 2 +1 2 = 73 

The unit eigenvectors are 


2 


1 


-1 

714 


76 


73 

1 

Unit Xi ,= 

J_ 

Unit Xi ,= 

1 

7l4 


76 


73 

3 


2 


1 

_7l4 


76 


LvsJ 


We observe that for the first unit eigenvector the sum of the squares is unity, that is, 

and the same is true for the other two unit eigenvectors in (5.83). 


(5.83) 


(5.84) 


5.6 Circuit Analysis with State Variables 

In this section we will present two examples to illustrate how the state variable method is used in 
circuit analysis. 


Example 5.10 


For the circuit of Figure 5.7, the initial conditions are i L (0 ) = 0, and v c (0 ) = 0.5 V. Use the 
state variable method to compute i L (t) and v c (t). 


v s (t) = u 0 (t) 



v c (t) 


Figure 5.7. Circuit for Example 5.10 


5-22 Signals and Systems with MATLAB ® Computing and Simulink ® Modeling, Fourth Edition 

Copyright ® Orchard Publications 




























Relationship between State Equations and Laplace Transform 


A: [-1 -1; 3/4 0] 

B: [1 0]’ 

C:[0 1] 

D: [ 0 ] 

Initial conditions: [0 0] 

For the Transfer Fen block Function Block Parameters dialog box we have entered: 

Numerator coefficient: [ 1 ] 

Denominator coefficient: [1 1 1 ] 


Step 1 


Step 2 









x' = Ax+Bu 


w 

y = Cx+Du 


State-Space 


PP 

Scope 1 

Bus 

Creator 1 










1 

- w 

s^+s+1 


Transfer Fen 


P^ 

Scope 2 


Bus 

Creator 2 


Figure 5.14- Models to show the equivalence between relations (5.126) and (5.127) 

After the simulation command is executed, both Scope 1 and Scope 2 blocks display the input 
and output waveforms shown in Figure 5.15. 



Figure 5.15. Waveforms displayed by Scope 1 and Scope 2 blocks for the models in Figure 5.14 
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Chapter 6 The Impulse Response and Convolution 


At 

e = 


-0.5e 1 + 1.5e 3t -2e‘ + 2e 3t 

3 -t 3 ~ 3t , , -t „ . -3t 

-e --e 1.5e -0.5e 


The impulse response is obtained from (6.5), Page 6-1, that is, 


At. 


h(t)= x(t) = e bu 0 (t) 


then, 


h(t)= x(t) = 


—t —t _'H 

-0.5e +1.5e -2e +2e 

3 -t 3 - 3t . . -t -3t 

-e --e 1.5e -0.5e 



4 

0 


-2e 1 + 6e 3t 


u 0 (t) = 

3 -t 3 ~ 3t 
-e - -e 




2 2 


3 t u o(t) (6.10) 


In Example 5.10, Chapter 5, Page 5-22, we defined 


and 

Then, 


x, = 1 , 


X, = V r 


or 


h(t) = x 2 = v c (t) - 1.5e l -1.5e 3t 


h(t) = v c (t) = 1.5(e ‘-e 3t ) 


( 6 . 11 ) 


Of course, this answer is not the same as that of Example 5.10, because the inputs and initial con¬ 
ditions were defined differently. 


6.2 Even and Odd Functions of Time 

A function f(t) is an even function of time if the following relation holds. 


f(-t) = f(t) 


( 6 . 12 ) 


that is, if in an even function we replace t with -t, the function f(t) does not change. Thus, poly¬ 
nomials with even exponents only, and with or without constants, are even functions. For 
instance, the cosine function is an even function because it can be written as the power series 


1 t 2 , t 4 t 6 
cost = 1-+-+ .. 

2! 4! 6! 


Other examples of even functions are shown in Figure 6.4. 
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Solutions to End-of-Chapter Exercises 


4. 




i(t)*v 2 (t) = J v 1 (t)v 2 (t-T)di = J 4te 


2(t - T) dx = 4e _2t f l xe 2T di 
o J o 


From tables of integrals, 


ax 

Jxe ax dx = ^y(ax-l) 


and thus 


2x, 


h(t)*v 2 (t) = (4e 2t ) 6 (2 4 T ^ 


= e _2t [e 2t (2t - 1) - e°(—1)] 


Check: 


0 -2t _2t 

= e (2t— 1 + e ) = 2t+e - 1 


v 1 (t)*v 2 (t)o Vj(s) ■ v 2 (s), V! (s) = 4 /s“,V 2 (s) = l/(s + 2) 

V,(s)-V 2 (s) = -J— = -±- 
s (s + 2) s + 2s 


syms s t; ilaplace(4/(s A 3+2*s A 2)) 

ans = 

2 *t-l+exp(-2 *t) 


5. 


To use the convolution integral, we must first find the impulse response. It was found in Exer¬ 
cise 1 as 


and with the given values, 


h(t) = i(t) = (1/L)C (R/L)t u 0 (t) 


h(t) = e *u 0 (t) 
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Chapter 6 The Impulse Response and Convolution 


u 0 (t) 


R 

-AA/V 
1 Q 

1H 



When the input is the unit step u 0 (t), 


» ou 

i(t )| V; =u n(t) = J U 0 (t-T)h(T)dT 



in - “O'- 
A 


U„(-T) 


1 



t ,0 


= J (1) • e T dt = -e T | 0 = e T | t = (1 -e l )u 0 (t) 



6 . 


o- 

+ 1H 

v m (t) 


R§ V out(t) 
1 Q 


We will first compute the impulse response, that is, the output when the input is the delta func¬ 
tion, i.e., v in (t) = 8(t). Then, by KVL 

di T 


and with i L = x 


L dt iL = 8(t) 


1 • x + 1 ■ x = 8(t) 
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Solutions to End-of-Chapter Exercises 


or 


x = - x + S(t) 


By comparison with x = Ax + bu, we observe that A = -1 and b = 1. 

From (6.5) 

h(t) - e At bu 0 (t) = e 1 • 1 = e 1 

Now, we compute v out (t) when v in (t) = u 0 (t) - u 0 (t - 1) by convolving the impulse response 
h(t) with this input v in (t), that is, v out (t) = v in (t)*h(t). The remaining steps are as in Exam¬ 
ple 6.5 and are shown below. 



7. 


R 


' outv 


From Exercise 6, 


Then, for this circuit, 


v„ = 


'(1 -(1 -e _t ) = e _t ) 


0 - 

+ 

Vin(t) 

1) - 

—Wv- 

1 Q 

L§ 

-0 

+ 

l v out(t) 

? 1 H 

= V L = V 

in- v R 


-t 

- e 

0 < t < 1 


‘(e - 1) 

t > 1 



.0- e *(e — 1) = (1 -e)e ‘ = -1.732e 1 


0 < t < 1 
t > 1 


The plot for the time interval 0 < t < 5 is shown below. 
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Evaluation of the Coefficients 


. 2 
sin x 



2n 

Figure 7.5. Graphical proof of 1 (sinmt) 


dt = re 


2 

COS X 



2 71 

Figure 7.6. Graphical proof of 1 (cosmt) 

J n 


dt - rt 


To evaluate any coefficient in (7.10), say b 2 , we multiply both sides of (7.10) by sin2t. Then, 

f(t)sin2t = ^a 0 sin2t + a, costsin2t + a 2 cos2tsin2t + a 3 cos3tsin2t + a 4 cos4tsin2t + ... 

2 

b 1 sintsin2t + b 2 (sin2t)“ + b 3 sin3tsin2t + b 4 sin4tsin2t + ... 

Next, we multiply both sides of the above expression by dt, and we integrate over the period 0 
to 2 jt . Then, 

271 i 271 2n 2it 

f(t)sin2tdt = -a 0 sin2tdt + ad costsin2tdt + a 2 cos2tsin2tdt 

~ n * n ^ n " n 


2n 

a 3 [ cos3tsin2tdt+ 

J o 


(7.11) 


2ti .271 2 n 

b,| sintsin2tdt + b 2 1 (sin2t)“dt + b, I sin3tsin2tdt + . 
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Chapter 7 Fourier Series 


that is, the shape of the negative half-cycle of the waveform is the same as that of the positive 
half-cycle, but inverted. 

We will test the most common waveforms for symmetry in Subsections 7.3.1 through 7.3.5 
below. 

7.3.1 Symmetry in Square Waveform 

For the waveform of Figure 7.7, the average value over one period T is zero, and therefore, 
a 0 = 0. It is also an odd function and has half-wave symmetry since —f(—t) = f(t) and 
-f (t + T/2) = f(t). 



An easy method to test for half-wave symmetry is to choose any half-period T/2 length on the 
time axis as shown in Figure 7.7, and observe the values of f(t) at the left and right points on the 
time axis, such as f(a) and f(b). If there is half-wave symmetry, these will always be equal but 
will have opposite signs as we slide the half-period T/2 length to the left or to the right on the 
time axis at non-zero values of f(t). 

7.3.2 Symmetry in Square Waveform with Ordinate Axis Shifted 

If in the square waveform of Figure 7.7 we shift the ordinate axis n/2 radians to the right, as 
shown in Figure 7.8, we will observe that the square waveform now becomes an even function, 
and has half-wave symmetry since f(-t) = f(t) and -f (t + T/2) = f(t). Also, a 0 = 0. 

Obviously, if the ordinate axis is shifted by any other value other than an odd multiple of n/2 , 
the waveform will have neither odd nor even symmetry. 
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The Exponential Form of the Fourier Series 


The waveform of Figure 7.30 is a rudimentary presentation of the capacitor voltage for the circuit 
of Figure 7.27. However, it will improve if we add a sufficient number of harmonics in (7.87). 

We can obtain a more accurate waveform for the capacitor voltage of Figure 7.27 with the Sim- 
ulink model of Figure 7.31. 



Pulse Transfer Fen Bus 

Generator Creator 


Figure 7.31. Simulink model for the circuit of Figure 7.27 
The input and output waveforms are shown in Figure 7.32. 



Figure 7.32. Input and output waveforms for the model of Figure 7.31 


7.8 The Exponential Form of the Fourier Series 

The Fourier series are often expressed in exponential form. The advantage of the exponential 
form is that we only need to perform one integration rather than two, one for the a n , and 

another for the b n coefficients in the trigonometric form of the series. Moreover, in most cases 
the integration is simpler. 

The exponential form is derived from the trigonometric form by substitution of 

jtot -joot 

coscot = --- (7.88) 

j oat — j cot 

sineot = ---— (7.89) 

J2 
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Chapter 8 The Fourier Transform 

8.2 Special Forms of the Fourier Transform 

The time function f(t) is, in general, complex, and thus we can express it as the sum of its real 
and imaginary parts, that is, as 

f(t) = f Re (t)+jf Im (t) (8.6) 

The subscripts Re and Im will be used often to denote the real and imaginary parts respectively. 
These notations have the same meaning as Re {f(t)} and Im {f(t)}. 

By substitution of (8.6) into the Fourier integral of (8.1), we obtain 

oo 00 

F(w) = J f Re (t)e~ J “ t dt + jj f Im (t)e~ JC0t dt (8.7) 

—00 —00 

and by Euler’s identity 


(8.9) 

( 8 . 10 ) 


00 00 

F(eo) = | [ f Re (t)coscot + f Im (t)sincot]dt - j| [ f Re (t)sincot-f Im (t)coscot]dt (8.8) 

—00 —00 

From (8.8), we see that the real and imaginary parts of F(eo) are 

oo 

F Re(®) = J [fReW COSCOt+ f Im (t)sin<Bt]dt 

—00 

and 

00 

F imO) = “J [ f Re (t)sincot-f Im (t)cos(ot]dt 

—00 

We can derive similar forms for the Inverse Fourier transform as follows: 

Substitution of (8.2) into (8.3) yields 

1 00 

f(t) = ^{ [ F Re ((°) + j F Im (®)] e JCOt d© 

and by Euler’s identity, 

Y oo 

f(t) = — [F Re (ro)cosrot-F Im (ro)sineot]deo 
2n J -oo 

Y oo 

+ j — f [ F Re (co)sincot + F Ira (©)coseot]deo 

^ —00 

Therefore, the real and imaginary parts of f(t) in terms of the Inverse Fourier transform are 


( 8 . 11 ) 


( 8 . 12 ) 


j oo 

f R e(t) = [ F Re (®) cos ®t-F Im (co) sina)t]dco 


(8.13) 


and 
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Special Forms of the Fourier Transform 


TABLE 8.3 Time Domain and Frequency Domain Correspondence (Refer also to Tables 8.4 - 8.7) 


f(t) 

F(©) 

Real 

Imaginary 

Complex 

Even 

Odd 

Real 



✓ 



Real and Even 

✓ 



✓ 


Real and Odd 






Imaginary 






Imaginary and Even 






Imaginary and Odd 







00 


F Re (©) = | f Re (t)cOS©tdt = 0 

—00 

fReM = odd 


and 


00 


F im(®) = “ 2 f f Re (t)sin©tdt 

J o 

fRe(t) = odd 


Therefore, if f Re (t) = odd, F(eo) is imaginary. 


( 8 . 20 ) 

( 8 . 21 ) 


To determine whether F(eo) is even or odd when f Re (t) = odd, we perform a test for evenness 
or oddness with respect to eo. Thus, substitution of -© for © in (8.21), yields 


oo oo 

F im(-®) = ~ 2 \ f Re (t)sin(-©)tdt = 2 f f Re (t)sin©tdt = -F Im (a) (8.22) 

J o J o 

Conclusion: If f(t) is real and odd, F(©) is imaginary and odd. We indicate this result in 
Table 8.4. 


TABLE 8.4 Time Domain and Frequency Domain Correspondence (Refer also to Tables 8.5 - 8.7) 


f(t) 

F(®) 

Real 

Imaginary 

Complex 

Even 

Odd 

Real 



✓ 



Real and Even 

✓ 



✓ 


Real and Odd 


✓ 



✓ 

Imaginary 






Imaginary and Even 






Imaginary and Odd 
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Chapter 8 The Fourier Transform 


Example 8.5 



1 . 

<=> 

2 J 


1 2 



syms t v w x; ft=t*exp(-t A 2); Fw=fourier (ft) 

Fw = 

-l/2*i*pi' s (1/2) *w*exp (- l/4*w' v 2 ) 

pretty(Fw) 


1/2 2 
- 1/2 i pi w exp(- 1/4 w ) 


Example 8.6 

-e~u 0 (t) + 38(t) «•-— + 3 

j©+ 1 

syms t v w x; fourier(sym('-exp(-t)*Heaviside(t)+3*Dirac(t)')) 

ans = 

-1/(l+i*w)+3 


(8.98) 


(8.99) 


Example 8.7 

u n (t) <=> tcS(co) + — 

J® 

syms t v w x; uO=sym('Heaviside(t)'); Fw=fourier(uO) 

Fw = 

pi*Dirac(w)-i/w 

We have summarized the most common Fourier transform pairs in Table 8.9. 


8.8 The System Function and Applications to Circuit Analysis 


We recall from Chapter 6 that, by definition, the convolution integral is 


We let 


00 

h(t )*u(t) = | u(t-T)h(t)dt 

—00 


g(t) = f(t)*h(t) 


( 8 . 100 ) 


( 8 . 101 ) 

( 8 . 102 ) 
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Chapter 8 The Fourier Transform 


We could also have used numerical integration with the integral 

f 


n • 2 

* sm x 


dx 


0 X 


thereby avoiding the integration by parts procedure. Shown below are the function 
fourierxfm2 which is saved as an .m file and the program execution using this function. 

function y2=fourierxfm2(x) 

x=x+(x==0)*eps; 

y2=(sin(x)./x). A 2; 

and after this file is saved, we execute the statement below observing that the limits of integra¬ 
tion are from 0 to n. 

value2=quad('fourierxfm2',0,pi) 

value2 = 

1.4182 


fplot('abs(2.*exp(-j.*w)*(sin(w)/w)) , ,[0 4*pi 0 2]) 
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Chapter 9 Discrete-Time Systems and the Z Transform 


and for m = 2, reduces to 


2{f[n + 2]} = z 2 F(z)-f [0]z 2 -f [l]z 


9.2.5 Multiplication by a" in the Discrete-Time Domain 


a f [n] <=> F - 
a 


Proof: 


~ 1 00 

{a n f [n]} = ^a n f[n]z- n = £ — n f[n]z' n = £ f [n] 


k = o 


k = 0 


k = 0 


9.2.6 Multiplication by e " aT in the Discrete-Time Domain 


e naT f [n] <=> F(e aT z) 


Proof: 


2 {e naT f [n]} = £ e 1,1,1 f [n]z 11 = £ f [n](e ai z) " = F(e ai z) 
k = 0 k = 0 

9.2.7 Multiplication by n and n 2 in the Discrete-Time Domain 


nf [n] <=> —z—F(z) 
dz 

2 d 2 d 2 

n“f [n] <=> z—F(z) + z —-F(z) 


dz 


dz“ 


Proof: 

By definition, 


F (z) = £ f[n]z~ 


n = 0 


and taking the first derivative of both sides with respect to z, we obtain 


A 00 . . 00 

—F(z) = £ (-n)f [n]z n = -z £ nf [n]z 11 


n = 0 


n = 0 


(9.16) 


(9.17) 


if z d 

n / z) 

- 

= F - 

v a J 

v a J 


(9.18) 


(9.19) 
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Properties and Theorems of the Z Transform 


Y(z) = 2 \ ^ x[m]h[n-m] l - ^ 


^ x[m]h[n-m] 


m = 0 J n = 0 L m = o 

and interchanging the order of the summation, we obtain 

00 00 00 00 

Y(z) = II x[m]h[n-m]z n = ^ x[m] h[n-m]z 

m = 0 L n = 0 -I m = 0 n = 0 

Next, we let k = n - m, then, n = k + m, and thus, 

00 00 /i \ 00 00 1 

Y(z) = x[m] h[k]z +m = ^ x[m]z m ^ h[k]z 


m = 0 n = 0 


m = 0 


n = 0 


or 


Y(z) - X(z) • H(z) 


(9.27) 


9.2.10 Convolution in the Discrete-Frequency Domain 

If f^n] and f 2 [n] are two sequences with 2 transforms Fj(z) and F 2 (z) respectively, then, 


fi[n]-f 2 [n] xF 1 (v)F 2 ^v ‘dv 


(9.28) 


where v is a dummy variable, and j) is a closed contour inside the overlap convergence regions 
for Xj(v) and X 2 (z/v). The proof requires contour integration; it will not be provided here. 


9.2.11 Initial Value Theorem 


Proof: 

For all n > 1 , as z —> go 


f [0] = limX(z) 


(9.29) 


~ n 1 „ n 

z =- > 0 


and under these conditions f [n]z n — > 0 also. Taking the limit as z — > oo in 

00 

F(z) = f[n]z n 


n = 0 
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The Inverse Z Transform 


num=[0 1 0 0]; % The coefficients of the numerator 

den=[1 -9/4 13/8 -3/8]; % The coefficients of the denominator 

fprintfC \n'); 

[num,den]=residue(num,den); % Verify the residues in (9.79) 
fprintf('r1 = %4.2f \t', num(1)); fprintf('p1 = %4.2f \t\ den(1));... 
fprintf('r2 = %4.2f \t\ num(2)); fprintf('p2 = %4.2f \t', den(2));... 
fprintf('r3 = %4.2f \t', num(3)); fprintf('p3 = %4.2f \t\ den(3)) 

rl = 8.00 pi = 1.00 r2 = -9.00 p2 = 0.75 r3 = 2.00 p3 = 0.50 

syms n z 

fn=2*(0.5) A n-9*(0.75) A n+8; % This is the answer in (9.80) 

Fz=ztrans(fn,n,z); simple(Fz) % Verify answer by first taking Z transform of f[n] 

ans = 

8 * z A 3/ (2*z-l)/(4*z-3)/ (z-1) 

iztrans(Fz) % Now, verify that Inverse of F(z) gives back f[n] 

ans = 

2*(1/2) A n-9*(3/4) A n+8 

We can use Microsoft Excel to obtain and plot the values of f [n] . The spreadsheet of Figure 9.6 
shows the first 25 values of n but only part of the spreadsheet is shown. 


n f[n] 

0.000 1.0000 
1.000 2.2500 
2.000 3.438 
3.000 4.453 
4.000 5.277 
5.000 5.927 
6.000 6.429 
7.000 6.814 
8.000 7.107 
9.000 7.328 
10.000 7.495 
11.000 7.621 
12.000 7.715 
13.000 7.786 
14.000 7.84 

15.000 7.88 

Figure 9.6. The discrete-time sequence f[n] = 2(0.5) n - 9(0.75) n + 8 for Example 9.4 


Discrete Time Sequence f[n] = 2(0.5)" - 9(0.75)" + 8 
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The Inverse Z Transform 


or 

([„] = i 6+ (1«/64)(0.75)° = 16 _i|3 (0 75) » 

(~0.25)(0.75) 2 9 

for n > 2. 

We can express f [n] for all n > 0 as 

f[n] - f5[n] + |5[n-l] + 16-^(0.75) n (9.96) 

where the coefficients of 8[n] and S[n- 1] are the residues that were found in (9.92) and (9.93) 
for n = 0 and n = 1 respectively at z = 0. The coefficient 28/9 is multiplied by 8[n] to 
emphasize that this value exists only for n = 0 and coefficient 4/3 is multiplied by 8[n- 1] to 
emphasize that this value exists only for n = 1 . 

Check with MATLAB: 

symsz n; Fz=(z A 3+2*z A 2+1)/(z*(z-1)*(z-0.75)); iztrans(Fz) 

ans = 

4/3 *charfcn[1] (n)+28/9*charfcn[0] (n)+16-163/9*(3/4) A n 
We evaluate and plot f [n] for the first 20 values. This is shown on the spreadsheet of Figure 9.9. 


n 

f[n] 


0 

1.000 


1 

3.750 


2 

5.813 

16 

3 

8.359 


4 

10.270 

12 

5 

11.702 

6 

12.777 


7 

13.582 

8 

8 

14.187 


9 

14.640 

4 

10 

14.980 


11 

15.235 


12 

15.426 

0 

13 

15.570 



Discrete Time Sequence of f[0] = 1, f[1] = 3.75, 
and f[n] = 16-(163/9)*(3/4) n for n >2 


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 


Figure 9.9. The discrete-time sequence for Example 9.7 


Example 9.8 

Use the inversion integral method to find the Inverse Z transform of 


Signals and Systems with MATLAB ® Computing and Simulink ® Modeling, Fourth Edition 9-35 
Copyright ® Orchard Publications 

































Chapter 10 The DFT and the FFT Algorithm 


F(z) = £f[n]z“ 


( 10 . 1 ) 


n = 0 


Its value on the unit circle of the z -plane, is 

JtoT 


F(e J ) - ^ f [n]e 


-njo)T 


( 10 . 2 ) 


n = 0 



This represents an infinite summation; it must be truncated before it can be computed. Let this 
truncated version be represented by 


(10.3) 


where N represents the number of points that are equally spaced in the interval 0 to 2n on the 
unit circle of the z -plane, and 

for m = 0, 1, 2, ..., N - 1. We refer to relation (10.3) as the N - point DFT of X [m] . 

The Inverse DFT is defined as 


(10.4) 


for n = 0, 1,2, . 

In general, the discrete frequency transform X [m] is complex, and thus we can express it as 

X [m] = Re{X [m]} + Im{X [m]} (10.5) 

for m = 0, 1, 2, ..., N - 1. 

Since 



J 71 n 27tmn • • 2rcmn 

- cos-ism- 

N J N 


we can express (10.3) as 

X [m] = ]Tx[n]e = ^ x [n] N 


N — 1 -j2ji— N-l 0 N-l t 

n v—i r n 27 rmn \. r i 27 rmn 

cos-j \ x n srn- 


n = 0 


n = 0 


n = 0 


N 


( 10 . 6 ) 


(10.7) 
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Chapter 10 The DFT and the FFT Algorithm 



Figure 10.2. Plots of x[ n] and |X[m]| values for Example 10.3 


On the plot of |X[m]| in Figure 10.2, the first value X[0] = 46.70 represents the DC component. 
We observe that after the X[8] = 0.10 value, the magnitude of the frequency components for the 
range 9 < m < 15, are the mirror image of the components in the range 1 < m < 7 . This is not a 
coincidence; it is a fact that if x [n] is an N -point real discrete time function, only N/2 of the fre¬ 
quency components of |X[m]| are unique. 

Figure 10.3 shows typical discrete time and frequency magnitude waveforms, for a N = 16 -point 
DFT. 
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The Fast Fourier Transform (FFT) 

The sequences of (10.65) and (10.66) cannot be decomposed further. They justify the statement 
made earlier, that each computation produces a vector where each component of this vector, for 
n = 1, 2, 3, ..., 7 is obtained by one multiplication and one addition. This is often referred to as a 
butterfly operation. 

Substitution of (10.65) and (10.66) into (10.60), yields 

F eveJ m ] = F evenl[ m ] + W N 2m F e ven2[ m ] (10.67) 

Likewise, F odd [m] can be decomposed into DFTs of length 2; then, F[m] can be computed from 

F ( m ) - F even( m ) + W N m F odd (m) m = 0, 1, 2, .... 7 (10.68) 

for N = 8 . Of course, this procedure can be extended for any N that is divisible by 2. 

Figure 10.5 shows the signal flow graph of a decimation in time, in-place FFT algorithm for 
N = 8 , where the input is shuffled in accordance with the above procedure. The subscript N in 
W has been omitted for clarity. 


Column 0 (x[n]) Column 1 (N /4) Column 2 (N / 2) Column 3 (N) X[m] 
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Chapter 10 The DFT and the FFT Algorithm 


and it is used with the nodes where two solid lines terminate on them. The number inside the 
circles denote the power of W N , and the minus (-) sign below serves as a reminder that the 
bracketed term of the second equation involves a subtraction. Using this graph and the above 
equations, compute the frequency component X [3]. Verify that this is the same as in Example 
10.5. 



6. Plot the Fourier transform of the rectangular pulse shown below, using the MATLAB fft func¬ 
tion. Then, use the ifft function to verify that the inverse transformation produces the rectan¬ 
gular pulse. 



f(t) / 

\ 

1 






-1 0 

1 


7. Plot the Fourier transform of the triangular pulse shown below using the MATFAB fft func¬ 
tion. Then, use the ifft function to verify that the inverse transformation produces the rectan¬ 
gular pulse. 
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Basic Analog Filters 


11.2.3 RLC Band-Pass Filter 

The RLC network shown in Figure 11.8 is a basic analog band-pass filter. We will derive expres¬ 
sions for its magnitude and phase. 

jcoL 1/jcoC 

+—'wnr*—1(—| + 

An R < V out 


Figure 11.8. Basic band-pass RLC network 
Application of the voltage division expression yields 

V = -—-• V 

out jcoL + 1/jcoC + R 111 

and thus 


G(jco) = XsHt = -R- = -- = - -MW - (ii.v) 

V in jcoL + 1/jcoC + R -co 2 LC + 1 + jeoRC co 2 - jco(R/L) - 1/LC 

There is no need to express relation (11.7) in magnitude and phase form. We will make use of 
the abs(x) and angle(x) MATLAB functions for the magnitude and phase plots. 

We will use the MATLAB script below to plot |G(jco)| versus radian frequency co. This is shown 
in Figure 11.9 where, for convenience, we let R = L = C = 1, and thus (11.7) simplifies to: 


G(jco) = ,~ j( ° (11.8) 

co" — j CO - 1 

w=0:0.02:100; magGjw=abs(-j.*w./(w. A 2-j.*w-1)); semilogx(w,magGjw);... 
xlabel('Frequency in rad/sec - log scale 1 ); ylabel('Magnitude of Vout/Vin 1 );... 
title('Magnitude Characteristics of basic RLC band-pass filter 1 ); grid 

The plot for the magnitude of (11.8) is shown in Figure 11.9. 

To plot the phase of (11.8), we use the MATLAB script below. 

w=0:0.02:100; phaseGjw=angle(-j.*w./(w. A 2-j.*w-1)).*180./pi; semilogx(w,phaseGjw);... 
xlabel('Frequency in rad/sec - log scale 1 ); ylabel('Phase of Vout/Vin - degrees 1 );... 
title('Phase Characteristics of basic RLC band-pass filter 1 ); grid 

The plot for the phase of (11.8) is shown in Figure 11.10. 


* This is the same network as (a) in Exercise 7, Chapter 4, Page 4-22. 
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Chapter 11 Analog and Digital Filters 

xlabel('Frequency in rad/sec - log scale 1 ); ylabel('Phase of Vout/Vin - degrees 1 );... 
title('Phase Characteristics of basic RLC band-elimination filter 1 ); grid 

The plot for the phase of (11.10) is shown in Figure 11.13. 


Phase Characteristics of basic RLC band-elimination filter 



Figure 11.13. Phase characteristics of the basic RLC band-elimination filter with R - L = C = 1 

Other passive filter networks are presented as end-of-chapter exercises in Chapter 4. Examples 
of active low-pass filters are presented in Figure 4.20, Chapter 4, Page 4-15, and at the end of this 
chapter. In the next section we introduce filter prototypes and use these as a basis for filter 
design. 

11.3 Low-Pass Analog Filter Prototypes 

In this section, we will use the analog low-pass filter as a basis. We will see later that, using trans¬ 
formations, we can derive high-pass and the other types of filters from a basic low-pass filter. We 
will discuss the Butterworth, Chebyshev Type I, Chebyshev Type II also known as Inverted Chebyshev, 
and Cauer also known as elliptic filters. 

The first step in the design of an analog low-pass filter is to derive a suitable magnitude-squared 
function A“(©), and from it derive a G(s) function such that 

A 2 (g>) = G(s)-G(-s)| s = jw (11.11) 

Since (j©)* = (-j©), the square of the magnitude of a complex number can be expressed as that 
number and its complex conjugate. Thus, if the magnitude is A, then 

A 2 (©) = |G(j©)| = G(j©)G*(j©) = G(j©) • G(-j©) (11.12) 

Now, G(j©) can be considered as G(s) evaluated at s = j©, and thus (11.11) is justified. Also, 
since A is understood to represent the magnitude, it needs not be enclosed in vertical lines. 
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Chapter 11 Analog and Digital Filters 


2nd Order Butterworth Low-Pass Filter Response 



Figure 11.20. Plot for the VC VS low-pass filter of Example 11.6 


We have used the MATLAB buttap function earlier to aid us in the design of Butterworth filters 
with the cutoff frequency normalized to 1 rad/s. We can also use the bode function to display 
both the (asymptotic) magnitude and phase plots. The following script will produce the Bode 
magnitude and phase plots for a two-pole Butterworth low-pass filter. 

[z,p,k]= buttap(2); % Specify a two-pole filter;... 

[b,a]=zp2tf(z,p,k); % Display in polynomial rational form;... 

w=0:0.01:4; [mag,phase]=bode(b,a,w);... 

b,a % Display b and a coefficients 

b = 

0 0 1 

a = 

1.0000 1.4142 1.0000 

num=[0 0 1];den=[1 sqrt(2) 1];... 

bode(num,den); title('Butterworth 2nd Order Low-Pass Filter 1 ); grid 

The Bode plots are shown in Figure 11.21. The frequency is displayed in rad/sec and the cutoff 
frequency normalized to 1 rad/s. 

We can also display the Bode plots with the frequency specified in Hz . This can be done with the 
MATLAB script below. 

h=bodeplot(tf(num,den));... 
setoptions(h,'FreqUnits', 'Hz'); grid 

The Bode plots with the frequency specified in Hz are shown in Figure 11.22. 
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Low-Pass Analog Filter Prototypes 


Type 2 Chebyshev Low-Pass Filter, k=3, 3 dB ripple in stop band 



11.3.4 Elliptic Analog Low-Pass Filter Design 

The elliptic, also known as Cauer filters, are characterized by the low-pass magnitude-squared 
function 

A 2 (co) = --i- (11.69) 

1 + R k (eo/eo c ) 

where R k (x) represents a rational elliptic function used with elliptic integrals. Elliptic filters have 
ripple in both the pass-band and the stop-band as shown in Figure 11.30. 


5-pole Elliptic Low Pass Filter 



Figure 11.30. Characteristics of an elliptic low-pass filter 
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Digital Filters 


Four signals separated by four digital filters 



Figure 11 .43. Frequency separations for the signals of Example 11.21 


In the following example, we will demonstrate the MATLAB filter function that is being used to 
remove unwanted frequency components from a function. But before we use the filter function, 
we must design a filter that is capable of removing those unwanted components. 


Example 11.22 

In Chapter 7, Subsection 7.4.4, Page 7-20, we found that the half-wave rectifier with no symme¬ 
try can be represented by the trigonometric Fourier series 


, A , A . A 

f(t) = — + — sint - — 

K 2 TT 


cos2t cos4t cos6t 


L 3 


15 


35 


cos8t 

63 


In this example, we want to filter out just the first 2 terms, in other words, to remove all cosine 
terms. To simplify this expression, we let A = 3 7t and we truncate it by eliminating all cosine 
terms except the cos2t term. Then, 

g(t) = 3 + 1.5sint-cos2t (11.99) 

The problem now reduces to design a low-pass digital filter, and use the filter command to 
remove the cosine term in (11.99). 

Solution: 

We will use a 6 - pole digital low-pass Butterworth filter because we must have a sharp transi¬ 
tion between the 1 and 2 rad/s frequency range. Also, since the highest frequency component 
is 2 rad/s, to avoid aliasing, we must specify a sampling frequency of co s = 4 rad/s. Thus, the 
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Digital Filters 


The analog and digital inputs and outputs are shown in Figures 11.45 and 11.46 respectively. 



Continuous Time t 


Figure 11.45. Continuous time input and output waveforms for Example 1122 



Figure 11.46. Discrete time input and output waveforms for Example 11.22 


We conclude this section with one more example to illustrate the use of the MATLAB find 
function. This function displays the subscripts of an array of numbers where a relational expres¬ 
sion is true. For example, 

x=-2:5; % Display the integers in the range -2 <= x < =5 

x = 

-2 -1 0 1 2 3 4 5 

k=find(x>0.8); % Find the subscripts of the numbers for which x > 0.8 
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Chapter 11 Analog and Digital Filters 


To implement this filter, we factor the numerator and denominator polynomials as 

H(z) = 0.5(1 + 0-6z 1 )(1 - 0.6Z 1 ) * 

(1 + 0.9z _1 )(l - 0.8z -1 ) 


(11.106) 


The model is shown in Figure 11.54, and the input and output waveforms are shown in Figure 
11.55. 


Series Form Realization - Second-Order Digital Filter - Simulation Time: 200 sec - Equ (11.106) 



Figure 11.54■ Model for Example 11.25 



Figure 11.55. Input and output waveforms for the model of Figure 11.54 


A demo model using fixed-point Simulink blocks can be displayed by typing 

fxpdemo_series_cascade_form 


* The combination of the of factors in parentheses is immaterial. For instance, we can group the factors as 

(1±m4) and (±z^£h oras (1±m4) and (A^h 

(l+0.9z ) (l-0.8z ) (l-0.8z ) (l+0.9z ) 
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Chapter 11 Analog and Digital Filters 


2. The op amp circuit below is a VCVS second-order band-pass filter whose transfer function is 

G( S ) = V °" t(s) = KFBWIs 

V in( S ) S 2 + [BW]s + (Oo 



Let eo 0 = center frequency, eo 2 = upper cutoff frequency, op = lower cutoff frequency, 
Bandwidth BW = eo 2 -©!, and Quality Factor Q = eo 0 /BW 

We can calculate the values of C ls C 2 , R l5 R 2 , R 3 , and R 4 to achieve the desired centered fre¬ 
quency cd 0 and bandwidth BW. For this circuit, 


R. = -2Q- 

CjOOoK 


R 2 = 


2Q 


R 3 = 


Ci(% 

1 


1 + 


7(K-1) 2 + 8Q 2 


1 + 1 


C 4 ©o Rl R 


R 4 — R 5 — 2R 3 

Using these relations, compute the appropriate values of the resistors to achieve center fre¬ 
quency f 0 = 1 KHz, Gain K = 10, and Q = 10. 

Choose the capacitors as Cj = C 2 = 0.1 pF. Plot |G(s)| versus frequency. 

Solution using MATLAB is highly recommended. 
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Appendix E Window Functions 




Figure E.2. Resultant amplitude response when the impulse response is abruptly truncated 

Therefore, we are seeking a suitable window function that will result in an amplitude response 
with a lower ripple in the stop band such as the one shown in Figure E.3. 


A(f) 



Figure E.3. Acceptable amplitude response 

Although an impulse response may converge uniformly, when it is multiplied by a rectangular 
window function it results in an undesirable amplitude response. Flowever, the rectangular win¬ 
dow function is the basis in studying window functions, so we begin with the description of the 
rectangular window function. 

E.2.1 Rectangular Window Function 

The rectangular window function is defined as 


f(t)rect = 1 for |t| < | 
- 0 otherwise 

and its time-domain waveform is as shown in Figure E.4. 



f(t) 






-t/2 t/2 


Figure E.4. Rectangular window function in time-domain 


* We recall that multiplication in the time domain corresponds to convolution in the frequency domain. 
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Appendix E Window Functions 



Figure E.13. Typical amplitude response for the Hanning window function 

Figure E.14 shows the normalized frequency domain plot for the Hanning window function cre¬ 
ated with the hann MATLAB function, and the script below. 

n=50; w=hann(n); [W,f]=freqz(w/sum(w),1,512,2); plot(f,20*log10(abs(W)));grid 



Figure E. 14■ Normalized frequency domain plot for the Hanning window function created with MATLAB 

We can also use the MATLAB Window Visualization Tool. The script below generates the plot 
shown in Figure E.15. 

L=50; wvtool(hann(L)) 
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Other Window Functions 

smaller side lobe levels. 

MATLAB provides the Kaiser window function as 

w=kaiser(L, beta) 

and this returns an L-point Kaiser window in the column vector w. 

Figure E. 24 shows the Kaiser functions for L = 50 with [3 = 1, [3 = 4, and P = 9. These func¬ 
tions were plotted with the MATLAB script below. 

plot([kaiser(50, 1), kaiser(50,4), kaiser(50,9)]); 



Figure E. 24- Kaiser window functions with L=50 and P = 1, P = 4 , and P = 9 

Figure E.25 shows the corresponding frequency domain plots for the Kaiser window functions in 
Figure E.24. These plots were created with the MATLAB script below. 

n=50; 

wl =kaiser(n, 1); w2=kaiser(n,4); w3=kaiser(n,9); 

[W1 ,f]=freqz(w1/sum(w1),1,512,2); [W2,f]=freqz(w2/sum(w2),1,512,2); 
[W3,f]=freqz(w3/sum(w3),1,512,2); plot(f,20*log10(abs([W1 W2 W3])));grid; 

E.3 Other Window Functions 

Table E.l lists the window functions we’ve discussed in the previous sections, as well as others 
along with the MATLAB function names and are included in the MATLAB Signal Processing 
Toolbox. 


* 


When the expression under the radical in (E.J2) is negative, the function can be expressed in terms of the 
hyperbolic sine function. 
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Appendix E Window Functions 


M 

Hi(z) = Y c m z m (E.16) 

m = -M 

Relation (E. 16) represents an FIR transfer function but it is non-causal because it includes posi¬ 
tive powers of z, and this implies that the filter would produce an output advanced in time with 
respect to the input, and this is, of course, impossible. We can be overcome this problem by intro¬ 
ducing a delay of M samples. Accordingly, we define the transfer function 

M -M 

H(z) = z- M Hl (z) = z- M Y c m z m = Y c m z " (E.17) 

m = -M m = M 

In (E.17) we let m-M = -i. Then, with the substitution c m _ i = a ; we obtain 


H(z) 


2M 


Z 


2M 


= z 


a i z 


(E-18) 


Let w m denote the coefficients of a particular window function and let c' m denote the coeffi¬ 
cients of the modified function, i.e., the given function multiplied by some window function. The 
modified coefficients are then computed from 


c = w c 
m mm 


(E.19) 


Example E. 1 

A low-pass FIR digital filter is to be designed using the Fourier series method. The desired ampli¬ 
tude response is 


A(f) = 1 for 0<f< 125 Hz 

0 elsewhere in the range 0 < f < f 0 


(E.20) 


The sampling frequency is 1 KHz and the impulse response is to be limited to 20 delays. Derive 
the transfer function using the Fourier series method. 

Solution: 

To make this low-pass filter an even function of frequency, we represent it as shown in Figure 
E.26 (a) in terms of the actual frequency. Let us express it in terms of the normalized frequency 
scale. We do this by considering the sampling frequency which is given as f s = 1 KHz. Accord¬ 
ingly, the folding frequency is f 0 = f s /2 = 1000/2 = 500 Hz, and in terms of the normalized 


* We recall that the folding frequency is the highest frequency that can be processed by a given discrete—time sys¬ 
tem with sampling frequency f s 
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Fourier Series Method for Approximating an FIR Amplitude Response 




Figure E.31. Normalized frequency plots for the triang window created with the MATLAB function firl 


Magnitude Response (dB) 



Figure E.32. Magnitude response for the triang window created with the MATLAB function fvtool 
b. Hanning 

We recall from (E.6) that the Hanning window function is defined as 

f(t) H ann = COS^Tlt/l) = + COS^ for |t| < | 

- 0 otherwise 
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Symbols 

% (percent) symbol in MATLAB A-2 

A 

abs(z) in MATLAB A-23 
active analog filter - see filter 
adjoint of a matrix - see matrix 
admittance 
capacitive 4-2 
inductive 4-2 
complex input 4-11 

algebraic constrain block in Simulink B-18 
aliasing 10-14 
all-pass filter - see filter 
all-pole approximation 11-21 
all-pole low-pass filter 
see filter - low-pass 
alternate form of the trigonometric 
Fourier series - see Fourier series 
alternate method of partial 

fraction expansion - see partial 
fraction expansion 
angle(z) MATLAB function A-23 
argument 11-2 
attenuation rate 11-12 
autoscale icon in Simulink B-12 
axis MATLAB command A-16 

B 

band-elimination filter - see filter 
band-elimination filter design 
see filter design 
band-limited signal 10-13 
band-pass filter - see filter 
band-pass filter design 
see filter design 
band-stop filter - see filter 
Bessel filter - see filter 
bilinear MATLAB function 11 -59 
bilinear transformation - see 

transformation methods for mapping 
analog prototype filters to digital filters 
bode MATLAB function 11-24 
box MATLAB command A-12 
buttap MATLAB function 11-17 
buttefly operation 10-23 
Butterworth analog low-pass 
filter design - see filter design 

C 

c2d MATLAB function 9-46 
capacitive admittance - see admittance 
capacitive impedance - see impedance 
cascade form realization - see digital filter 
Category I FFT algorithm - see 
FFT algorithm 

Category II FFT algorithm - see 
FFT algorithm 


Cauchy’s residue theorem 
see residue theorem 
Cauer filter - see elliptic filter 
Cayley-Flamilton theorem 5-11 
characteristic equation 5-19 
cheblap MATLAB function 11-35 
cheb2ap MATLAB function 11-38 
Chebyshev filters - see filter 
Chebyshev Type I analog low-pass 
filter design - see filter design 
Chebyshev Type I filters - see filter 
Chebyshev Type I low-pass filter 
magnitude-square function 11-26 
- see filter 

Chebyshev Type II analog low-pass 
filter design - see filter design 
Chebyshev Type II filters - see filter 
circuit analysis with Laplace transforms 4-1 
circuit analysis with state variables 5-22 
circuit transformation from time 
to complex frequency 4-1 
clc MATLAB command A-2 
clear MATLAB command A-2 
cofactor of a matrix - see matrix 
collect(s) MATLAB symbolic function 3-12 
column vector in MATLAB A-19 
command screen in MATLAB A-1 
Command Window in MATLAB A-1 
commas in MATLAB A-8 
comment line in MATLAB A-2 
Commonly Used Blocks in Simulink B-7 
complex admittance - see admittance 
complex conjugate in MATLAB A-4 
complex conjugate pairs 3-5 
complex impedance - see impedance 
complex number C-2 
complex numbers in MATLAB A-3 
complex poles 3-5 
Complex to Magnitude-Angle 
block in Simulink C-7 
computation of the state 
transition matrix 5-11 
computation of the Z Transform with 
contour integration - see Z transform 
Configuration Parameters 
in Simulink B-12 

congugate of a matrix - see matrix 
conj(A) MATLAB function D-9 
conjugate of a complex number C-3 
conjugate time and frequency functions 
of the Fourier transform - see 
Fourier transform - properties of 
constant function - Fourier transform of 
see Fourier transform - properties of 
Contents Pane in Simulink B-7 
contour integral 9-20 
conv MATLAB function A-7 
convolution in the complex frequency 
domain - see Laplace transform 
properties of 

convolution in the discrete-frequency 
domain - see Z transform - properties of 


convolution in the discrete-time 

domain - see Z transform - properties of 
convolution in the time domain property 
of the Laplace transform - see 
Laplace transform - properties of 
convolution integral defined 6-8 
graphical evaluation of 6-8 
convolution property of the Fourier 
transform - see Fourier transform 
properties of 
Cooley and Tukey 10-18 
cosine function - Fourier transform of 
see Fourier transform of 
common functions 

cosco 0 t u 0 (t) function - Fourier transform of 
see Fourier transform of 
common functions 
Cramer’s rule D-17 

D 

d2c MATLAB function 9-47 
data points in MATLAB A-14 
decade - definition of 11-12 
decimation in frequency - see FFT algorithm 
decimation in time - see FFT algorithm 
deconv MATLAB function A-6 
default color in MATLAB A-15 
default line in MATLAB A-15 
default marker in MATLAB A-15 
default values in MATLAB A-12 
delta (impulse) function 
definition of 1-11 
doublet 1-14 

Fourier transform of - see Fourier 
transform of common functions 
higher order 1-14 
nth-order 1-14 
sampling property of 1-12 
sifting property of 1 -13 
triplet 1-14 

demo in MATLAB A-2 
DeMoivre's theorem 11-15 
derivation of the Fourier transform from 
the Laplace Transform 8-25 
determinant of a matrix - see matrix 
determinant of order 2 - see matrix 
DFT - common properties of 
even time function 10-9 
even frequency function 10-9 
frequency convolution 10-13 
frequency shift 10-12 
linearity 10-10 
odd time function 10-9 
odd frequency function 10-9 
time convolution 10-12 
time shift 10-11 
DFT - definition of 10-1 
N-point 10-2, 10-16 

diagonal elements of a matrix - see matrix 
diagonal matrix - see matrix 
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differentiation in complex frequency 
domain property of the Laplace 
transform - see Laplace transform 
- properties of 

differentiation in time domain property 
of the Laplace transform - see 
Laplace transform - properties of 
differentiation property of the Fourier 
transform 8-12 - see Fourier 
transform - properties of 
digital filter 11-1, 11-51, 11-70 

Finite Impulse Response (FIR) 11-52 
FIR 11-52 
HR 11-51 

Infinite Impulse Response (HR) 11-51 
realization of 

Direct Form I 11-70 
Direct Form II 11-71 
cascade (series) form 11 -73 
non-recursive 11-52 
parallel form 11-70 
recursive 11-52 
series (cascade) form 11 -73 
Digital Filter Design Simulink block 11 -78 
digital filter design with Simulink 11 -70 
dimpulse MATLAB function 9-29 
Dirac MATLAB function 1-20 
Direct Form I realization - see digital filter 
realization of 

Direct Form II realization - see digital filter 
realization of 

direct term in MATLAB 3-4 
discontinuous function - definition of 1-2 
Discrete Fourier Transform (DFT) 10-1 
discrete impulse response 9-40 
discrete-time system transfer function 9-40 
discrete unit step function 9-3 
discrete-time exponential sequence 9-16 
discrete-time systems 9-1 
discrete-time unit ramp function 9-18 
discrete-time unit step function 9-14 
Display block in Simulink B-18 
display formats in MATLAB A-31 
distinct eigenvalues - see eigenvalues 
distinct poles - see poles 
division of complex numbers C-4 
dot operator in MATLAB 
division with A-21 
exponentiation with A-21 
multiplication with A-20 
double-memory technique 
see FFT Algorithm 
doublet - see delta function 

E 

Editor Window in MATLAB A-1 
Editor/Debugger in MATLAB A-1 
eig(x) MATLAB function 5-17 
eigenvalues 
distinct 15-1 

multiple (repeated) 5-15 
eigenvector 5-19 

e"'“* u 0 (t) Fourier transform of - see Fourier 
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transform of common functions 
element-by-element operation in MATLAB 
division A-21 
exponentiation A-21 
multiplication A-20 
elements of the matrix - see matrix 
ellip MATLAB function 11-40 
elliptic filter - see filter 
elliptic filter design - see filter design 
eps in MATLAB A-22 
Euler’s identities C-5 
even functions 6-4, 7-33 
even symmetry - see Fourier 
series - symmetry 
Excel's Analysis ToolPak 10-5 
exit MATLAB command A-2 
expand MATLAB symbolic function 3-10 
exponential form of complex numbers C-5 
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see Fourier series 
exponential order function - 
definition of 2-2 
eye(n) MATLAB function D-7 
eye(size(A)) MATLAB function D-7 
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factor(s) MATLAB symbolic function 3-4 
Fast Fourier Transform (FFT) 10-1, 10-17 
FDA Tool Digital Filter Design 
Simulink block 11-82 
FFT algorithm 
Category I 10-19 
Category II 10-20 
decimation in frequency 10-20 
decimation in time 10-20 
double-memory technique 10-20 
in-place 10-20 
natural input-output 10-20 
FFT definition of 10-1, 10-17 
fft(x) MATLAB function 10-5, 11-68 
Figure Window in MATLAB A-13 
filter - see also digital filter 
active 

high-pass 4-22, 4-32 
low-pass 4-22, 4-32, 4-23, 4-35 
all-pass 11-1, 11-94 
all-pole 11-21 

band-elimination 11-1, 11-8 
band-pass 11-1, 11-7 
band-stop - see band-elimination 
Bessel 11-95 
Chebyshev 11-10 
Inverted 11-38 

magnitude-square function 11-26 
prototype 11-10 
Type 111-25 
Type II 11 -38 
elliptic 11-39 

high-pass 4-22, 4-30, 11-1,11-4 
low-pass 4-22, 4-30, 11-1,11-2 
low-pass analog filter 
prototypes 11-10 
maximally flat 11-14 (footnote) 


notch (band-elimination) 11-8 
phase shift 11-1. 11-94 
RC high-pass 11-4 
RC low-pass 11 -2 

RLC band-elimination 4-22, 4-31,11-8 
RLC band-pass 4-22, 4-31,11-7 
filter design - see also digital filter 
band-elimination 11-41 
band-pass 11-41 

Butterworth analog low-pass 11-14 
Chebyshev 
Type 111-25 
Type II 11-38 
elliptic 11-39 
high-pass 11-41 
low-pass 11-14 
filter MATLAB function 11-63 
final value theorem in Z transform 
see Z transform - properties of 
final value theorem in Laplace transform 
see Laplace transform - properties of 
find MATLAB function 11-68 
Finite Impulse Response (FIR) digital filter 
see digital filter 
FIR - see digital filter 
first harmonic - see Fourier series 
harmonics of 
first-order circuit 5-1 
first-order simultaneous 
differential equations 5-1 
Flip Block command in Simulink B-11 
format in MATLAB A-31 
fourier MATLAB command 8-33 
Fourier series 

exponential form of 7-31 
method used in window functions E-17 
trigonometric form of 7-2, 7-10 
alternate form of 7-25 
Fourier series coefficients - evaluation of 
numerical evaluation using Excel 7-46 
numerical evaluation using MATLAB® 7-47 
Fourier series of common waveforms 
full-wave rectifier 7-20, 7-24 
half-wave rectifier 7-17, 7-20 
square waveform 

with even symmetry 7-9, 7-13 
with odd symmetry 7-8, 7-12 
sawtooth 7-9, 7-15 
triangular 7-9, 7-16 
Fourier series - harmonics of 
first 7-1,7-10 
second 7-1,7-10 
third 7-1,7-10 
Fourier series - symmetry 
even 7-7 

half-wave 7-7, 7-34 
odd 7-7 

quarter-wave 7-7 (footnote) 
types of 7-7 

Fourier integral - see Fourier transform 
Fourier transform 
definition of 8-1 
inverse of 8-1 
special forms of 8-2 



Fourier transform - properties of 
area under f(t) 8-15 
area under F(co) 8-15 
conjugate time and frequency 
functions 8-13 
constant function 8-18 
frequency convolution 8-15 
frequency differentiation 8-13 
frequency shifting 8-11 
imaginary time functions - Fourier 
transform of 8-6 
linearity 8-9 

Parseval's theorem 8-16 
real time functions - Fourier 
transform of 8-3 
symmetry 8-9 
time convolution 8-14 
time differentiation 8-12 
time integration 8-13 
time scaling 8-10 
time shifting 8-11 
Fourier Transform derivation from 
Laplace transform 8-25 
Fourier transform of common functions 
cosro 0 t 8-19 
cosro 0 t u 0 (t) 8-24 
delta (8(t) and S(t-a)) 8-18 
e''“ ot 8-19 
e >ot u 0 (t) 8-24 
signum (sgn(t)) 8-20 
sinco 0 t 8-20 
sinco 0 t u 0 (t) 8-25 
unit step (u 0 (t)) 8-22 

Fourier transform of common waveforms 
combined rectangular pulses 8-29 
cosine within a rectangular pulse 8-30 
shifted rectangular pulse 8-28 
symmetrical rectangular pulse 8-27 
periodic time functions 8-31,8-32 
fourth-order all-pole low-pass filter 
see filter, low-pass 
fplot in MATLAB A-27 
frequency convolution in DFT 
see DFT - common properties of 
frequency convolution in Fourier 
transform of - see Fourier transform 

- properties of 

frequency differentiation in Fourier 
transform of - see Fourier transform 

- properties of 
frequency shift in DFT 

see DFT - common properties of 
frequency shift in Fourier transform 
see Fourier transform - properties of 
frequency shift in Laplace transform 
see Laplace transform - properties of 
freqz MATLAB function 11-57 
full rectification waveform - Laplace 
transform of - see Laplace transform 
of common waveforms 
full-wave rectifier - Fourier series of - see 
Fourier series of common waveforms 
Function Block Parameters in Simulink B- 
function files in MATLAB A-26 
fundamental frequency 7-1 


fzero MATLAB function A-26 

G 

Gain block in Simulink B-9, B-18 
gamma function 2-15 
Gaussian elimination method D-19 
generalized factorial function 2-15 
Gibbs phenomenon 7-24 
grid MATLAB command A-12 
gtext MATLAB command A-13 

H 

half-wave rectifier - Fourier series of - see 
Fourier series of common waveforms 
half-wave symmetry 

see Fourier series - symmetry 
Heavyside MATLAB function 1-14 
help in MATLAB A-2 
Hermitian matrix - see matrix 
higher order delta functions - see delta 
function 

high-pass filter - see filter 

high-pass filter design - see filter design 

I 

identity matrix - see matrix 
ifft(x) MATLAB function 10-5 
ifourier MATLAB function 8-33 
HR - see digital filter 
ilaplace MATLAB function 3-4 
imag(z) MATLAB function A-23 
imaginary axis - definition of C-2 
imaginary number - definition of C-2 
imaginary time functions 8-6 

see Fourier transform - properties of 
impedance 
capacitive 4-2 
inductive 4-2 
complex input 4-8, 4-9 
improper integral - definition of 2-15 
improper rational function - 
definition of 3-1,3-13 
impulse function - see delta function 
impulse invariant method - see 

transformation methods for mapping 
analog prototype filters to 
digital filters 

increments between points 
in MATLAB A-14 

inductive admittance - see admittance 
inductive impedance - see impedance 
infinite impulse response - see digital filter 
initial value theorem in Z transform 
see Z transform - properties of 
initial value theorem in Laplace transform 
see Laplace transform - properties of 
in-place FFT algorithm 10-20 
see FFT algorithm 

integration in frequency in Laplace transform 
10 see Laplace transform - properties of 
integration in time in Laplace transform 
see Laplace transform - properties of 


Inverse Fourier transform 
see Fourier transform 
Inverse Laplace transform 2-1 
see Laplace transform 
inverse of a matrix - see matrix 
Inverse Z transform - see Z transform 
inversion integral 9-32 
Inverted Chebyshev filter - see filter 


j operator C-1 

L 

L’ Hopital’s rule 2-16 
laplace MATLAB function 2-27 
Laplace integral - see Laplace transform 
Laplace transform 
definition of 2-1 
Inverse of 2-1,3-1 
Laplace transform - properties of 
convolution in the complex 
frequency domain 2-12 
convolution in the time domain 2-11 
differentiation in complex 
frequency domain 2-6 
differentiation in time domain 2-4 
final value theorem 2-10 
frequency shift 2-4 
initial value theorem 2-9 
integration in 

complex frequency domain 2-8 
time domain 2-6 
linearity 2-3 
scaling 2-4 
time periodicity 2-8 
time shift 2-3 

Laplace transform of common waveforms 
full-rectified 2-32 
half-rectified 2-26 
linear segment 2-23 
pulse 2-23 

rectangular periodic 2-25 
sawtooth 2-32 
triangular 2-24 

Laplace transform of common functions 
transform of e' at coscot u 0 (t) 2-22 
transform of e' at sincot u 0 (t) 2-21 
transform of e' at u 0 (t) 2-19 
transform of coscot u 0 (t) 2-20 
transform of S(t) 2-18 
transform of 5(t-a) 2-18 
transform of sincot u 0 (t) 2-20 
transform of t" u 0 (t) function 2-15 
transform of t" e' al u 0 (t) 2-19 
transform of u 0 (t) 2-14 
transform of u,(t) 2-14 
leakage 10-13 

left shift in discrete-time domain 
see Z transform - properties of 
Leibnitz’s rule 6 

lims= MATLAB command A-27 

line spectra 7-36 

linear difference equation 9-38 


IN-3 



linearity in DFT 

see DFT - common properties of 
linearity in discrete-time domain 
see Z transform - properties of 
linearity property in Fourier transform 
see Fourier transform - properties of 
linearity property in Laplace transform 
see Laplace transform - properties of 
linspace MATLAB command A-14 
In (natural log) A-13 
log in MATLAB A-13 
log(x) MATLAB function A-13 
Iog10(x) MATLAB function A-13 
log2(x) MATLAB function A-13 
loglog MATLAB command A-13 
long division of polynomials 9-36 
lower triangular matrix - see matrix 
low-pass analog filter prototypes - see filter 
low-pass filter - see filter 
Ip2bp MATLAB function 11-43 
Ip2bs MATLAB function 11-43 
Ip2hp MATLAB function 11-43 
Ip2lp MATLAB function 11-43 

M 

magnitude-squared function 11-10 
main diagonal of a matrix - see matrix 
Math Operations in Simulink B-10 
MATLAB Demos A-2 
MATLAB's Editor/Debugger A-1 
matrix (matrices) 
adjoint of D-21 
cofactor of D-13 
conformable for addition D-2 
conformable for subtraction D-2 
conformable for multiplication D-4 
congugate of D-8 
definition of D-1 
determinant D-10 
minor of D-13 
non-singular D-21 
singular D-21 
diagonal D-2, D-6 
diagonal elements of D-2 
elements of D-1 
Hermitian D-9 
identity D-7 
inverse of D-22 
left division in MATLAB D-25 
multiplication in MATLAB A-18 
power series of 5-9 
scalar D-7 
size of D-7 
skew-Hermitian D-9 
skew-symmetric D-9 
square D-1 
symmetric D-8 
trace of D-2 
transpose of D-8 
triangular 
lower D-6 
upper D-6 
zero D-2 

matrix left division in MATLAB - see matrix 


matrix multiplication in MATLAB - see matrix 
matrix power series - see matrix 
maximally flat filter - see filter 
mesh(x,y,z) MATLAB function A-17 
meshgrid(x,y) MATLAB command A-17 
m-file in MATLAB A-2, A-26 
minor of determinant - see matrix 
MINVERSE Excel function D-27 
MMULT Excel function D-27 
modulated signals 8-12 
multiple eigenvalues - see eigenvalues 
multiple poles - see poles 
multiplication by an in discrete-time domain 
see Z transform - properties of 
multiplication by e' naT in discrete-time 
domain - see Z transform - properties of 
multiplication by n in discrete-time domain 
see Z transform - properties of 
multiplication by n 2 indiscrete-time domain 
see Z transform - properties of 
multiplication of complex numbers C-3 

N 

NaN in MATLAB A-26 
natural input-output FFT algorithm 
see FFT algorithm 
network transformation 
resistive 4-1 
capacitive 4-1 
inductive 4-1 

non-recursive realization digital filter 
see digital filter 

non-singular determinant - see matrix 
normalized cutoff frequency 11-15 
notch filter - see filter 
N-point DFT - see DFT - definition of 
nth-order delta function - see delta function 
numerical evaluation of Fourier coefficients 
see Fourier series coefficients 
Nyquist frequency 10-13 

O 

octave defined 11-12 
odd functions 6-5, 7-34 
odd symmetry - see Fourier 
series - symmetry 
orthogonal functions 7-2 
orthogonal vectors 5-19 
orthonormal basis 5-19 

P 

parallel form realization - see digital filter 
Parseval’s theorem - see 

Fourier transform - properties of 
partial fraction expansion 3-1,3-2, 9-25 
alternate method of 3-15 
method of clearing the fractions 3-15 
phase angle 11-2 
phase shift filter - see filter 
picket-fence effect 10-14 
plot MATLAB command A-10 
polar form of complex numbers C-6 


polar plot in MATLAB A-24 
polar(theta.r) MATLAB function A-23 
poles 3-1 
complex 3-5 
distinct 3-2 

multiple (repeated) 3-8 
poly MATLAB function A-4 
polyder MATLAB function A-7 
polynomial construction from 
known roots in MATLAB A-4 
polyval MATLAB function A-6 
pre-sampling filter 10-13 
pre-warping 11-55 
proper rational function - 
definition of 3-1, 11-11 
properties of the DFT 

see DFT - common properties of 
properties of the Fourier Transform 
see Fourier transform - properties of 
properties of the Laplace Transform 
see Laplace transform - properties of 
properties of the Z Transform 
see Z transform - properties of 

Q 

quarter-wave symmetry 7-7 (footnote) 
quit MATLAB command A-2 

R 

radius of absolute convergence 9-3 
ramp function 1-9 
randn MATLAB function 11-68 
Random Source Simulink block 11 -80 
rationalization of the quotient C-4 
RC high-pass filter - see filter 
RC low-pass filter - see filter 
real axis C-2 
real number C-2 
real(z) MATLAB function A-23 
rectangular form C-5 
rectangular pulse expressed in terms 
of the unit step function 1-4 
recursive realization digital filter 
see digital filter 
region of 

convergence 9-3 
divergence 9-3 

relationship between state equations 
and Laplace Transform 5-30 
residue 3-2, 9-41 

residue MATLAB function 3-3, 3-12 
residue theorem 9-20, 9-21 
right shift in the discrete-time domain 
see Z transform - properties of 
RLC band-elimination filter - see filter 
RLC band-pass filter - see filter 
roots of polynomials in MATLAB A-3 
roots(p) MATLAB function 3-6, A-3 
round(n) MATLAB function A-24 
row vector in MATLAB A-3 
Runge-Kutta method 5-1 
running Simulink B-7 
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s 

sampling property of the delta function 
see delta function 
sampling theorem 10-13 
sawtooth waveform - see Laplace 
transform of common waveforms 
sawtooth waveform - Fourier series of 
see Fourier series of 
common waveforms 
scalar matrix - see matrix 
scaling property of the Laplace transform 
see Laplace transform - properties of 
Scope block in Simulink B-12 
script file in MATLAB A-2, A-26 
second harmonic - see Fourier series 
harmonics of 

semicolons in MATLAB A-8 
semilogx MATLAB command A-12 
semilogy MATLAB command A-12 
series form realization - see digital filter 
Shannon’s sampling theorem 
see sampling theorem 
shift of f[n] u 0 [n] in discrete-time domain 
see Z transform - properties of 
sifting property of the delta function 
see delta function 
signal flow graph 10-23 
signals described in math form 1-1 
signum function - see Fourier transform 
of common functions 
simout To Workspace block 
in Simulink B-12 

simple MATLAB symbolic function 3-7 
Simulation drop menu in Simulink B-12 
simulation start icon in Simulink B-12 
Simulink icon B-7 
Simulink Library Browser B-8 
sine function - Fourier transform of 
see Fourier transform of 
common functions 
singular determinant - see matrix 
Sinks library in Simulink B-18 
sino) 0 t u 0 (t) Fourier transform of - see 
Fourier transform of common functions 
size of a matrix - see matrix 
skew-Hermitian matrix - see matrix 
skew-symmetric matrix - see matrix 
special forms of the Fourier transform 
see Fourier transform 
spectrum analyzer 7-36 
square matrix - see matrix 
square waveform with even symmetry - see 
Fourier series of common waveforms 
square waveform with odd symmetry - see 
Fourier series of common waveforms 
ss2tf MATLAB function 5-33 
stability 11-13 

start simulation in Simulink B-12 
state equations 

for continuous-time systems 5-1 
for discrete-time systems 9-45 
state transition matrix 5-9 
state variables 

for continuous-time systems 5-1 


for discrete-time systems 9-45 
State-Space block in Simulink B-12 
state-space equations 

for continuous-time systems 5-1 
for discrete-time systems 9-45 
step function - see unit step function 
step invariant method - see trans¬ 
formation methods for mapping analog 
prototype filters to digital filters 
stop-band filter - see filter 
string in MATLAB A-16 
subplots in MATLAB A-18 
summation in the discrete-time Domain 
see Z transform - properties of 
symmetric matrix - see matrix 
symmetric rectangular pulse expressed 
as sum of unit step functions 1-6 
symmetric triangular waveform expressed 
as sum of unit step functions 1-6 
symmetry - see Fourier series - symmetry 
symmetry property of the Fourier transform 
see Fourier transform - properties of 
system function - definition of 8-35 

T 

Taylor series 5-1 
text MATLAB command A-14 
tf2ss MATLAB function 5-33 
theorems of the DFT 10-10 
theorems of the Fourier Transform 8-9 
theorems of the Laplace transform 2-2 
theorems of the Z Transform 9-3 
third harmonic - see Fourier 
series - harmonics of 
time convolution in DFT 

see DFT - common properties of 
time integration property of the Fourier 
transform - see Fourier 
transform - properties of 
time periodicity property of the Laplace 
transform 2-8 - see Laplace 
transform - properties of 
time scaling property of the Fourier 
transform - see Fourier 
transform - properties of 
time shift in DFT 

see DFT - common properties of 
time shift property of the Fourier transform 
see Fourier transform - properties of 
time shift property of the Laplace transform 
see Laplace transform - properties of 
title(‘string’) MATLAB command A-12 
trace of a matrix - see matrix 
Transfer Fen block in Simulink 4-17 
Transfer Fen Direct Form II 
Simulink block 11-71 
transfer function of 

continuous-time systems 4-13 
discrete-time systems 9-38 
transformation between 
s and z domains 9-22 
transformation methods for mapping 
analog prototype filters to digital filters 
Impulse Invariant Method 11-52 


Step Invariant Method 11-52 
Bilinear transformation 11-53 
transpose of a matrix - see matrix 
Tree Pane in Simulink B-7 
triangular waveform expressed in terms 
of the unit step function 1-4 
triplet - see delta function 
Tukey - see Cooley and Tukey 

U 

unit eigenvectors 5-19 
unit impulse function (8(t)) 1-8 
unit ramp function (u^t)) 1-8 
unit step function (u 0 (t)) 1-2 
upper triangular matrix - see matrix 
using MATLAB for finding the Laplace 
transforms of time functions 2-27 
using MATLAB for finding the Fourier 
transforms of time function 8-33 

V 

Vandermonde matrix 10-18 
Vector Scope Simulink block 11 -84 

W 

warping 11-54 
window functions 
Blackman E-12 

Fourier series method for approximating 
an FIR amplitude response E-17 
Hamming E-9, E-31 
Hanning E-7, E-27 
Kaiser E-14, E-35 

other used as MATLAB functions E-15 
rectangular E-2 
triangular E-5, E-23 

Window Visualization Tool in MATLAB E-4 

X 

xlabel MATLAB command A-12 

Y 

ylabel MATLAB command A-12 

Z 

Z transform 

computation of with contour 
integration 9-20 
definition of 9-1 
Inverse of 9-1,9-25 
Z transform - properties of 
convolution in the discrete 
frequency domain 9-9 
convolution in the discrete 
time domain 9-8 
final value theorem 9-10 
initial value theorem 9-9 
left shift 9-5 
linearity 9-3 
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multiplication by a" 9-6 
multiplication by e' naT 9-6 
multiplication by n 9-6 
multiplication by n 2 9-6 
right shift 9-4 
shift of f[n] u 0 [n] 9-3 
summation 9-7 

Z Transform of discrete-time functions 
cosine function cosnaT 9-16 
exponential sequence e' naT u 0 [n] 
9-16, 9-21 

geometric sequence a" 9-11 
sine function sinnaT 9-16 
unit ramp function nu 0 [n] 9-18,9-21 
unit step function u 0 [n] 9-14, 9-20 
zero matrix - see matrix 
zeros 3-1,3-2 

zp2tf MATLAB function 11-17 
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